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Abstract 

In a recently published paper Ciufolini reports on the so far per- 
formed tests aimed at the detection of the general relativistic gravito- 
magnetic Lense-Thirring effect in the gravitational field of the Earth by 
means of the analysis of the laser-ranged data of the existing LAGEOS 
and LAGEOS II geodetic satellites. In this paper we will critically dis- 
cuss his claims by showing that the total error, mainly due to the 
systematic bias due to the mismodelling in the static and time- varying 
parts of the multipolar expansion of the Newtonian terrestrial gravi- 
tational potential, is larger than that claimed by Ciufolini. E.g., the 
systematic error due to the mismodelling in the static part of the geopo- 
tential in the tests performed with the EGM96 Earth gravity model and 
the combination involving the nodes of LAGEOS and LAGEOS II and 
the perigee of LAGEOS II realistically amounts to more than 80% (1- 
a): the claimed total uncertainty, including also the non-gravitational 
perturbations which especially affect the perigee of LAGEOS II, is, 
instead, 20-25%. The claimed accuracy in the more precise tests per- 
formed with the 2nd generation CHAMP-only EIGEN2 Earth gravity 
model and a combination involving the nodes of LAGEOS and LA- 
GEOS II over 10 years is 18%. With numerical simulations we will 
show that, instead, it is < 51% (l-c) if the impact of the secular vari- 
ations of the even zonal harmonics over a so long observational time 
span (~ 14%) is accounted for. 
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1 Introduction 

Recent years have seen increasing efforts aimed to directly detecting various 
phenomena connected to the general relativistic gravitomagnetic field ^ 
121 El Ej of the rotating Earth. It should be noted that, according to K. 
Nordtvedt [SJ, the multidecadal analysis of the Moon'orbit by means of the 
Lunar Laser Ranging (LLR) technique yields a comprehensive test of the 
various parts of order C(c -2 ) of the post-Newtonian equation of motion. The 
existence of gravitomagnetism as predicted by the Einstein's General Theory 
of Relativity would, then, be indirectly inferred from the high accuracy of 
the lunar orbital reconstruction. In jSJ the same arguments are applied to 
the radial motion of the LAGEOS satellite. 

The extraordinarily sophisticated and expensive Gravity Probe B (GP- 
B) mission |Hj has been launched in April 2004; it is aimed at the detection 
of the gravitomagnetic precession of the spins [S] of four superconducting 
gyroscopes carried onboard at a claimed accuracy of 1% or better. 

The Lense-Thirring effect on the orbital motion of a test particle ^U] 
could be measured by analyzing the orbital data of certain Earth artificial 
satellites with the Satellite Laser Ranging (SLR) technique Up to now, 
the only performed tests are due to Ciufolini and coworkers. 

In this paper we will analyze the latest results presented in |12j from 
a critical point of view in order to show that the claimed accuracies are 
optimistic. 

2 The Lense-Thirring effect on the orbit of a test 
particle and the strategy to measure it 

The gravitomagnetic field of a spinning mass of proper angular momentum 
J induces tiny secular precessions on the longitude of the ascending node f2 
and the argument of pericentre 1 a; of a test particle ^3 El HI EHj 

■ 2GJ . 6GJ cos i 

LT _ C 2 a 3(l- e 2)3/2' WLT - " c 2 a 3(! _ e 2)3/2' ( ) 

where G is the Newtonian constant of gravitation, c is the speed of light 
in vacuum, a, e and i are the semimajor axis, the eccentricity and the in- 
clination, respectively, of the test particle's orbit. In the terrestrial space 
environment the gravitomagnetic precessions are very small: for the geode- 
tic SLR LAGEOS satellites, whose orbital parameters are listed in Table 

1 In their original paper Lense and Thirring use the longitude of pericentre vj — + cu. 
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Table 1: Orbital parameters (a semimajor axis, e eccentricity, i inclination) 
of the existing LAGEOS and LAGEOS II satellites and their Lense-Thirring 
node precessions Clu? in m as yr -1 . 



Satellite 


a km) 


e 


i° 


$7 L t (mas yr" 1 ) 


LAGEOS 


12270 


0.0045 


110 


31 


LAGEOS II 


12163 


0.0135 


52.64 


31.5 



they amount to a few tens of milliarcseconds per year (mas yr -1 in the 
following) 

The extraction of the Lense-Thirring precessions from the orbit data 
analysis is very difficult due to a host of competing classical orbital pertur- 
bations of gravitational [13 HH1 H3 UH] and non-gravitational [T31 1201 I2T1 1221 
l23| l24*j origin which have various temporal signatures and are often quite 
larger than the relativistic signal of interest. The most insidious ones are the 
perturbations which have the same temporal signature of the Lense-Thirring 
precessions 2 , i.e. secular trends. Indeed, whatever the length of the adopted 
observational time span T b s is, they cannot be fitted and removed from the 
time series without removing the relativistic signal as well. Then, it is of 
the utmost importance to assess as more accurately and reliably as possible 
their aliasing impact on the measurement of the Lense-Thirring effect. 

It turns out that the perigees of the LAGEOS-like satellites are severely 
affected by the non-gravitational perturbations, contrary to the nodes. More- 
over, since the non-conservative forces depend on the structure, the shape 
and the rotational status of the satellite their correct modelling is not a 
trivial task and, as we will see later, introduces large uncertainties in the 
correct assessment of the error budget in some of the performed gravitomag- 
netic tests. 

2.1 The gravitational error 

The even (I = 2,4,6...) zonal (m = 0) harmonic coefficients Jt of the mul- 
tipolar expansion of the Earth's gravitational potential, called geopotential, 

2 Also the perturbations which grow quadratically in time are, of course, very danger- 
ous. Those induced by the secular variations of the even zonal harmonics of the Earth's 
geopotential fall in this category, as we will see in detail in Section E-i.2.11 Time-dependent 
periodic perturbations with periods longer than the observational time span may also be 
insidious because they would resemble superimposed linear trends |15| . 
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induce secular precessions on the node and the perigee of any near-Earth 
artificial satellite [251 which, of course, depend only on its orbital configura- 
tion and are independent of its physical structure. Such aliasing effects are 
many orders of magnitude larger than the Lense-Thirring precessions; the 
precision with which the even zonal harmonics are known in the currently 
available Earth gravity models [23113111113120112111311231131 would yield 
errors amounting to a significant fraction of the Lense-Thirring precessions 
or even larger. 

Even more dangerous are the perturbations induced by the secular vari- 
ations of the low degree even zonal harmonics Jg, I = 2, 4, 6 |35U36| . Indeed, 
such perturbations grow quadratically in time if the shifts in mas are con- 
sidered and linearly in time if the rates in mas yr _1 are considered. Their 
impact on the orbital elements of the LAGEOS satellites have been worked 
out in [37]. It turns out that, by using the results of [SHI; the errors induced 
by J2 would amount to 8%, 14% and 5.4% for the nodes of LAGEOS and 
LAGEOS II and the perigee of LAGEOS II, respectively, over an observa- 
tional time span T Q b s of just one year at 1— a level. This clearly shows that 
it would be impossible to analyze single orbital elements. 

The time-dependent periodic perturbations are less dangerous 

because if their periods are shorter than the adopted observational time span 
they can be fitted and removed from the time series. The most insidious 
tidal perturbation is that induced by the even zonal constituent which has 
a period of 18.6 years and whose nominal impact on the orbital elements 
of the LAGEOS satellites amounts to thousands of mas [15] . However, it 
turns out that it does not affect the observables which have been adopted for 
the performed Lense-Thirring tests because its main component is of degree 
1 = 2 and order m = 0. 

2.1.1 The linear combination approach 

The problem of reducing the impact of the mismodeling in the even zonal 
harmonics of the geopotential with the currently existing satellites can be 
coped in the following way [38] . 

Let us suppose we have at our disposal N (N> 1) time series of the resid- 
uals of those Keplerian orbital elements which are affected by the geopoten- 
tial with secular precessions, i.e. the node and the perigee: let them be ^> A , 
A=LAGEOS, LAGEOS II, etc. Let us write explicitly down the expressions 

3 Also the subtle non-gravitational Yarkovsky-Rubincam force, which is due to the 
interaction of the Earth's electromagnetic IR radiation with the physical structure of the 
LAGEOS satellites, induces secular effects on their nodes and perigees |21|. 
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of the observed residuals of the rates of those elements Stp^ in terms of 
the Lense-Thirring effect V>l T , °^ N-l mismodelled classical secular preces- 
sions ijj^SJe induced by those even zonal harmonics whose impact on the 
measurement of the gravitomagnetic effect is to be reduced and of the re- 
maining mismodelled phenomena A which affect the chosen orbital element 



fo/4 g = V^ T/ u LT + Yl ^ 6Je + AA ' A = LAGEOS, LAGEOS II,... 

N— 1 terms 

(2) 

The parameter 4 /xlt is equal to 1 in the General Theory of Relativity and 
in Newtonian mechanics. The coefficients tp^ are defined as 

and have been explicitly worked out for the node and the perigee up to 
degree i = 20 in they depend on some physical parameters of the 

central mass (GM and the mean equatorial radius R) and on the satellite's 
semimajor axis a, the eccentricity e and the inclination i. We can think 
about eq. (J2J) as an algebraic nonhomogeneuous linear system of N equa- 
tions in N unknowns which are //lt and the N-l 5Jf. solving it with respect 
to /Ult allows to obtain a linear combination of orbital residuals which is 
independent of the chosen N-l even zonal harmonics. In general, the or- 
bital elements employed are the nodes and the perigees and the even zonal 
harmonics cancelled are the first N-l low-degree ones. 

This approach is, in principle, very efficient in reducing the impact of 
the systematic error of gravitational origin because all the classical preces- 
sions induced by the static and time-dependent parts of the chosen N-l Ji 
do not affect the combination for the Lense-Thirring effect. Moreover, it 
is flexible because it can be applied to all satellites independently of their 
orbital configuration, contrary to the butterfly configuration in which the 
cancellation of the even zonal harmonics can be achieved only for supple- 
mentary orbital planes and identical orbital parameters. Apart from the first 
orbital element which enters the combination with 1, the other elements are 
weighted by multiplicative coefficients Cj(a, e, %) ^ 1 which are built up with 
ip.i and, then, depend on the orbital elements of the considered satellites. 
Their magnitude is very important with respect to the non-gravitational per- 
turbations, which in general are not cancelled out by the outlined method, 

4 It can be expressed in terms of the PPN 7 parameter as flur = (1 + t)/2. 
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and to the other time-dependent perturbations of gravitational origin with 
I ^ 2, 4, 6, ..,m 7^ 0. Values smaller than 1 for the Cj coefficients are, in 
general, preferable because they reduce the impact of such uncancelled per- 
turbations. It is important to note that the order with which the orbital 
elements enter the combination is important: indeed, while the systematic 
error due to the even zonal harmonics of the geopotential remains unchanged 
if the orbital elements of a combination are exchanged, the coefficients q do 
change and, consequently, also the non-gravitational error. The best results 
are obtained by choosing the highest altitude satellite as first one and by 
inserting the other satellites in order of decreasing altitudes. 

This method was explicitly adopted for the first time in [HE] with the 
nodes of the LAGEOS satellites and the perigee of LAGEOS II. The obtained 
combination is 

rrSLAGEOS , r , rrSLAGEOS II , „ r - LAGEOS II ,, fin 9 ( A\ 

where c\ = 0.304, C2 = —0.350 and 60.2 is the slope, in mas yr _1 , of the 
expected gravitomagnetic linear trend. Eq. Q is insensitive to the first two 
even zonal harmonics J2 and J4. It has been used in [lUj when the level of 
accuracy of the JGM3 26 and EGM96 [22j Earth gravity models, available 
at that time, made it necessary to consider a combination of observables 
which is independent of errors in both J2 and J4. 

In view of the great improvements in the Earth gravity field modelling 
with the CHAMP and, especially, GRACE missions an extensive 
search for alternative combinations has been subsequently performed |43 | I37 | 
EHEHIEHEZI- In [QUEUES] the following combination has been proposed 5 

^LAGEOS + fc^LAGEOS II „ (g) 

where k\ = 0.546 and 48.2 is the slope, in mas yr , of the expected gravit- 
omagnetic linear trend. It has been adopted for the tests performed in |12j 
with the 2nd generation CHAMP-only EIGEN2 Earth gravity model [50] 
and the 1st generation GRACE-only GGM01S [33] Earth gravity model. 
Eq. © allows to cancel out the first even zonal harmonic J2. 

The possibility of using only the nodes of the LAGEOS satellites in view of the im- 
provements in the Earth gravity models from GRACE has been propsed for the first time 
in |42| . although without quantitative details. In [12| it seems that Ciufolini refers to it 
as a proper own result with his reference [6] which includes |4U| of the present work and 
an announced paper; |40| is not concerned with eq. ||KJ| because it deals with eq. @ and 
its analysis by means of EGM96. Moreover, Iorio retains the e-mails in which he passed 
to Ciufolini, with whom he was long in contact, the combination of eq. @ along with the 
estimates of the systematic error obtained with EIGEN2 and is disposed to make them 
publicly available on request. 



7 



3 The performed Lense-Thirring tests with the 
LAGEOS satellites 

The only performed tests aimed at the detection of the Lense-Thirring pre- 
cessions of eq. (pQ) in the gravitational field of the Earth with the exist- 
ing LAGEOS satellites have been performed, up to now, by Ciufolini and 
coworkers. They have used the node-node-perigee combination of eq. (jlj) 
[101 IT2"] and the node-node combination of eq. ((SJ) ^2] ■ 

In |12j it is claimed that "[...] the Lense-Thirring effect exists and its 
experimental value, [...], fully agrees with the prediction of general relativity." 
in regard to both the tests with the EGM96 and EIGEN-2 Earth gravity 
models. In this Section we will disprove such statements. 

The main objections to the results presented in these works can be sum- 
marized as follows 

• Ciufolini has not performed tests by varying the length of the adopted 
observational time span, running backward and forward the initial 
epoch of the analysis, varying the secular rates of the even zonal har- 
monics in order to check their impact over different time spans, using 
different Earth gravity models in order to obtain a scatter plot of the 
obtained results. 

• The total error budget has been underestimated, especially the sys- 
tematic error of gravitational origin. E.g., the impact of the secular 
variations of the even zonal harmonics of the geopotential, which may 
become a very limiting factor over time spans many years long as those 
used, has not been addressed. Almost always 1— a results have been 
presented without any explicit indication of this fact. 

3.1 The node-node-perigee tests 

The combination of eq. Q) has been analyzed by using the EGM96 [22] 
Earth gravity model over 4 years in [IU] and over 7.3 years in |12j . The 
claimed total error budget amounts to 20-25% over 4 years and to 20% over 
7.3 years. 

3.1.1 The gravitational error 

The impact of the remaining uncancelled even zonal harmonics of the geopo- 
tential J&,Js, Jio, ... on eq. @ has been estimated by Ciufolini and cowork- 
ers with the full covariance matrix of EGM96 in a root-sum-square calcula- 
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tion. In [3U] and, six years later, in ^5] it is claimed to be < 13%. Apart 
from the fact that this is a 1 — a level estimate, in (35], as later acknowledged 
in a number of papers [3H1 EE HH EH 133 EH- the use of the full covariance 
matrix of EGM96 has been questioned. Indeed, it has been noted that in 
the EGM96 solution the recovered even zonal harmonics are strongly recip- 
rocally correlated; it seems, e.g., that the 13% value for the systematic error 
due to geopotential is due to a lucky correlation between Jq and J$ which 
are not cancelled by eq. (J3J). The point is that, according to j3H], nothing 
would assure that the covariance matrix of EGM96, which is based on a 
multi-year average that spans the 1970, 1980 and early 1990 decades, would 
reflect the true correlations between the even zonal harmonics during the 
particular time intervals of a few years adopted in the analyses by Ciufolini 
and coworkers. Then, a more conservative, although pessimistic, approach 
would be to consider the sum of the absolute values of the errors due to 
the single even zonal as representative of the systematic error induced by 
our uncertainty in the terrestrial gravitational field according to EGM96 
[33 OH- In this case we would get a conservative upper bound of 83% (1-0"). 
If a root-sum-square calculation is performed by neglecting the correlations 
between the even zonals a 45% 1-a error is obtained [161 1451 13T1 147j . 

3.1.2 The non-gravitational error 

Another important class of systematic errors is given by the non-gravitational 
perturbations which affect especially the perigee of LAGEOS II. The main 
problem is that it turned out that their interaction with the structure of 
LAGEOS II changes in time due to unpredictable modifications in the phys- 
ical properties of the LAGEOS II surface (orbital perturbations of radia- 
tive origin, e.g. the solar radiation pressure and the Earth albedo) and in 
the evolution of the spin dynamics of LAGEOS II (orbital perturbations of 
thermal origin induced by the interaction of the electromagnetic radiation 
of solar and terrestrial origin with the physical structure of the satellites, in 
particular with their corner-cube retroreflectors) . Moreover, such tiny but 
insidious effects were not entirely modelled in the GEODYN II software at 
the time of the analysis of [401 112j , so that it is not easy to correctly and 
reliably assess their impact on the total error budget of the measurement 
performed during that particular time span. According to the evaluations 
in [2j| , the systematic error due to the non-gravitational perturbations over 
a time span of 7 years amounts to almost 28%. However, according to [48] . 
their impact on the measurement of the Lense-Thirring effect with the nodes 
of LAGEOS and LAGEOS II and the perigee of LAGEOS II is, in general, 
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quite difficult to be reliably assessed. 

So, by adding quadratically the gravitational and non-gravitational er- 
rors of we obtain for the systematic uncertainty SfJr™ ematlc ^ 54% jf we 
assume a 45% error due to geopotential. The sum of the absolute values of 
the errors due to gepotential added quadratically with the non-gravitational 
perturbations would yield a total systematic error of 5 ^^ tematlc ^ 88%. It 
must be noted that the latter estimate is rather similar to those released 
in |48j . Note also that they are 1-a evaluations. Moreover, it should be 
considered that the perigee of LAGEOS II is also sensitive to the eclipses 
effect on certain non-gravitational perturbations. Such features are, gener- 
ally, not accounted for in all such estimates. An attempt can be found in 
|19| in which the impact of the eclipses on the effect of the direct solar radi- 
ation pressure on the LAGEOS-LAGEOS II Lense-Thirring measurement 
has been evaluated: it should amount to almost 10% over an observational 
time span of 4 years. 

3.2 The node-node tests 

In this Section we will deal with the node- node combination of eq. (J5J). 
Such observable only cancels out the gravitational bias of the first even 
zonal harmonic J2, but has the great advantage of discarding the perigee of 
LAGEOS II and its insidious non-gravitational perturbations. 

3.2.1 The gravitational error 

In |12| the node-node combination of eq. (JSJ) has been analyzed with the 2nd 
generation CHAMP-only EIGEN2 03 and the 1st generation GRACE-only 
GGM01S [33] Earth gravity models over a time span of almost 10 years. 

In the impact of the static part of the geopotential, according to the 
CHAMP-only EIGEN2 Earth gravity model, is evaluated as 18% (1-a root- 
sum-square covariance calculation), 22% (1-a root-sum-square calculation) 
and 37% (1-a upper bound). Ciufolini reports 18% obtained in a root-sum- 
square fashion with the full covariance matrix of EIGEN2 for which the same 
remarks as for EGM96 holds. Moreover, he does not consider the fact that 
EIGEN2 is only based on six months of data and that the released sigmas 
of the even zonal harmonics of low degree, which are the most relevant in 
this kind of analyses with the LAGEOS satellites, are rather optimistic, as 
explicitly pointed out in [3j]] and acknowledged in |37| . In regard to the 
GGM01S model, the covariance matrix was not publicly released. Ciufolini 
correctly presents a 19% which is the 1— a upper bound obtained in |37j . 
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However, GGM01S is only based on 111 days of data. 

In our opinion the author's conclusion "We conclude, using the Earth 
gravity model EIGEN-2S, that the Lense-Thirring effect exists and its ex- 
perimental value, fj, = 0.98 ±0.18, fully agrees with the prediction of general 
relativity" is optimisitc. Indeed, he claims that in his 18% total error bud- 
get all the error sources are included. Ciufolini neglects the impact of the 
time-dependent gravitational perturbations on eq. ©• Indeed, they may 
turn out to be a serious limiting factor mainly due to the secular variations 
of the even zonal harmonics 6 Jg. Indeed, eq. © allows to cancel out J2, 
but is sensitive to J4, Jq,-.-, as pointed out in [50] . The uncertainties in the 
J 1 are still quite large: see Table 1 of jHHl- From it the values of Table |2 can 
be inferred. 

Table 2: Weighted means and standard deviations from Table 1 of (SEl of 
the secular rates of the first three even zonal harmonics in units of 10 -11 





e = 2 


e = 4 


£ = 6 


Jl 


-2.113 


-0.6992 


-0.3594 


J 


0.0810 


0.2029 


0.1765 



On the other hand, their impact on the Lense-Thirring measurement 
grows linearly in time 7 . Indeed, the mismodelled shift, in mas, of eq. ® 
due to the secular variations of the uncancelled even zonal harmonics can 
be written as 

£ (r> L / GEOS + ^ L / GEOS n) a -frl s , (6) 

£=4 

where the coefficients are <9^ c iass/c^ and have explicitly been calculated 
up to degree £ = 20 in ^2^3- It must be divided by the gravitomagnetic 
shift, in mas, of eq. © over the same observational time span 

(O^geos + fc^AGEOS 11) ^ = mas yr -i Tobs _ (7) 

6 The problem of the secular variations of the even zonal harmonics in post-Newtonian 
tests of gravity with LAGEOS satellites has been quantitatively addressed for the first 
time in |51|. In regard to the Lense-Thirring measurement with eq. ©, it has been, 
perhaps, misunderstood in |37|. 

7 For a possible alternative combination which would cancel out the first three even 
zonal harmonics along with their temporal variations see |46l I47| . 
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By assuming, e.g., a j 4 = 0.6 x 1CT 11 yr _1 and aj g = 0.5 x 10~ n yr _1 [3"o] . 
it turns out that the percent error on the combination eq. (JSJ) grows linearly 
with T Q bs and would amount to 1% over one year at 1 — a level. This means 
that, over 10 years, their impact is ~ 10% (1-c). In Section 13.2.21 we will 
quantitatively support this evaluation. 

3.2.2 The impact of the secular variations of the even zonal har- 
monics: a quantitative estimate 

Here we describe a numerical experiment aimed at a quantitative evaluation 
of the impact of J^. 

The first step consists in simulating the time series of <5f2 LAGEOS + 
A;i50 LAGEOSn in order to obtain the qualitative and quantitative features of 
Figure 4 of ^21 • It refers to EIGEN2 and shows the raw residual time series 
with a straight line which fits it. The post-fit residuals amounts to 12 mas. 
In our model, called Input Model (IM), we include 

• LT= Surt with Slt = 48.2 mas yr _1 . Lense-Thirring trend as pre- 
dicted by the General Theory of Relativity in order to simulate the 
fact that the residuals of LAGEOS and LAGEOS II have been built 
up by dealing with the gravitomagnetic force as a totally unmodelled 
feature 

. ZONDOT^ £t 4 {r} (ft LAGEOS + Cl ft LAGEOS n ) (i) t 2 . Quadratic 

term due to the Ji according to Table |^1 The numbers {r} are ran- 
domly generated from a normal distribution with mean zero, variance 
one and standard deviation one. Note that EIGEN2 does not solve 
for J i. In ^2] there are no details about the values included in the 
dynamical force models of the orbital processor; thus we treat the sec- 
ular rates of the even zonal harmonics as unmodelled features fully 
absorbed by the residuals. 

• ZONALS= p (jHq) Surt. Linear trend with a slope of x% of the Lense- 
Thirring signal. For EIGEN2 x = 37 (sum of the absolute values of the 
individual errors) is assumed. The number p is randomly generated as 
for the {r} 

. TIDE= £{a c }cT4 c cos [(^) t + {/ c }] + ZM<JA S sin [(^) t + {/,}] . 
Set of various tidal perturbations of known periods P. For the impact 
of such kind of perturbations on the orbits of the LAGEOS satellites 
see |14] . The sets of numbers {a c }, {a s }, {/ c }, {f s } are randomly gen- 
erated as p and the {r} 
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• NOISE. White gaussian noise with variable amplitude which simulates 
the observational errors of the laser-ranged measurement 

The full IM used in our analysis is thus 



We include in our model the possibility of varying the length of the time 
series T b s , the temporal step At which simulates the orbital arc length, 
the amplitude of the noise and of the mismodelling in the perturbations 
and the initial phases of the sinusoidal terms in order to simulate different 
initial conditions and uncertainties in the dynamical force models of the 
orbital processors. More precisely, the magnitude of the mismodelling in 
the various effects is randomly varied within the currently accepted ranges 
(1-cr) by using random numbers generated from a normal distribution with 
mean zero, variance one and standard deviations one. The same also holds 
for ZONALS because it is impossible to know, a priori, the sign of the slope 
of the residual trend due to the even zonal harmonics. The so built IM 
represents the basis of our subsequent analyses. 

With the so built IM we perform a set of 5000 runs by randomly varying 
the initial phases, the noise and the mismodelling amplitudes within the 
accepted intervals in order to simulate a wide range of initial conditions and 
measurement errors which can occur in the real world. The length of the 
time series is keep fixed at T ^ s = 9.5 years. In every run we fit the IM with 
a straight line only (LF) and record the obtained slope \x with the related 
error 5(i. Then, we calculate the averages of [i and 5fj> and the standard 
deviation of the mean S. For EIGEN2, i.e. x = 37, we obtain 



Also the post-fit residuals are calculated. Figure ^ shows the complete IM, 
its straight line fit compared with the nominal Lense-Thirring trend and the 
post-fit residuals for a given set of randomly chosen initial conditions. The 
RMS post-fit amounts to 12.7 mas. The averaged RMS post-fit is 11 mas. 
This shows that our procedure represents a realistic starting point for our 
analyses. 

A first interesting result is that the departure of the measured slope 
fi from the nominal gravitomagnetic slope /Ult> which is included in IM, 
amounts to ~ 30% on average, while for Ciufolini is 2% only (/i = 0.98). 

In order to evaluate the impact of the secular rates of the even zonal 
harmonics in every run we also fit the IM with a quadratic polynomial (QF) 



IM = LT + ZONDOT + ZONALS + TIDE + NOISE. 



(8) 




(9) 
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Simulated time series 5£1 LAGE0S +k1 5£2 LAGEOS 




012345678 



Time (years) 



012345678 

Time (years) 



Figure 1: Simulated time series, straight line fit and post-fit residuals for 
T bs = 9.5 years and At = 15 days. The RMS of the post-fit residuals is 
12.7 mas. The slope of the trend simulating the impact of the mismodelled 
even zonal harmonics has been fixed to 37% of the Lense-Thirring effect, 
according to EIGEN2. 
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and compare the so obtained slope //qf with the slope obtained in LF /Ulf- 
Note that procedure is analogous to that adopted in ^H] for the periodic 
perturbations. On average, the difference between the two slopes, i.e. the 
systematic error due to J4 , Jq , amounts to 



A*LT 



14% 



(10) 



of the Lense-Thirring effect. This result holds for 1 — a. It is important 
to note that Figure ^ has been obtained by assuming that the combined 
residuals of the LAGEOS satellites absorb the quadratic signature of Ji 
according to Table [21 at l-cr, i.e. IM also includes ZONDOT. Nonetheless, it 
is difficult to discern the parabolic signal which, instead, is present and does 
affect the recovery of the slope. This means that a simple visual inspection 
of the plots of the combination of eq. (jSJ) cannot be considered conclusive 
about the effect of J4, Jq. 

Another important point is that the combination of eq. © cannot be 
used to reliably constrain the zonals' rates by measuring a jf. Indeed, it 
turns out that the errors 5Q in the quadratic parameters of QF are always 
larger than the estimated values themselves Q and their mean < SQ/Q > 
over a given set of 5000 runs amounts to ~ 260% with a standard deviation 
of the mean of 68%. Moreover, these figures change for different sets of 5000 
runs. 

If we repeat the same numerical experiments for T b s = 9.5 years with 
x = 18 (1-<t root-mean-square full covariance calculation) and x = 22 (l-cr 
root-mean-square variance calculation) the situation does not substantially 
change 



20%, 



'0. 



Ill) 



Then, a more conservative l-cr estimate of the total systematic error of 
the measurement of the Lense-Thirring effect with the combination of eq. 
(jSJ) and the EIGEN2 Earth gravity model is 



6fJ toM error < 5J% 



(12) 



While the forthcoming solutions from CHAMP and, especially, GRACE 
will be able to improve the static part of the terrestrial gravitational po- 
tential, i.e. the Je, it is not so for their secular rates Jg. This fact sets for 
the systematic error of gravitational origin a sort of threshold below which 
it will not be possible to go unless much more accurate determinations of 
J4, Jq will be available. 
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4 Conclusions 

In this paper we have performed a detailed critical analysis of the reliability 
and robustness of the so far performed tests aimed to the detection of the 
Lense-Thirring effect in the gravitational field of the Earth with the existing 
or proposed LAGEOS satellites. 

We can summarize our conclusions as follows 

• In regard to the node-node-perigee LAGEOS-LAGEOS II combina- 
tion, the claimed 20 — 25% total accuracy obtained with the EGM96 
Earth gravity model, still presented in is not realistic because of 
the impact of the non-gravitational perturbations on the perigee of 
LAGEOS II and the mismodelling in the even zonal harmonics of the 
geopotential whose 1 — a upper bound is 83%. 

• In regard to the node-node LAGEOS-LAGEOS II combination of eq. 
([5]). extensive numerical tests have been performed in order to quan- 
titatively assess the impact of the uncancelled secular variations of 
the even zonal harmonics on the proposed measurement of the Lense- 
Thirring effect. A simulated time series curve has been fitted with a 
straight line and a quadratic polynomial and the so obtained slopes 
have been compared. This procedure has been repeated over 5000 runs 
performed by randomly varying the initial phases, the noise and the 
mismodelling level within the currently accepted ranges of the simu- 
lated signal. It turns out that the bias due to J4 and Jq over 9.5 years 
amounts to ~ 14% on average. This yields an upper bound of the total 
systematic error of the performed tests with EIGEN2 of ~ 51%. 

• Alternative combinations involving the use of existing laser-ranged tar- 
gets other than the LAGEOS satellites should be analyzed. The most 
promising combination is, in principle, the one that involves the nodes 
of LAGEOS, LAGEOS II, Ajisai and Jason-1 [HI EH- It cancels out 
the first three even zonal harmonics J2,Ja,J& along with their tem- 
poral variations at the price of introducing the relatively huge non- 
gravitational perturbations on Jason-1 which, however, should have a 
time-dependent periodic signature with short periodicities. Accord- 
ing to the recently released combined CHAMP+GRACE+terrestrial 
gravimetry/altimetry EIGEN-CG01C Earth gravity model the 
systematic error due to the remaining even zonal harmonics would 
amount to 0.7 % (root-sum-square calculation) and 8 1.6% (upper bound) 

8 The 1-a errors for the node-node LAGEOS-LAGEOS II combination of eq. are 5% 



16 



at 1— a |46j . The forthcoming more accurate and robust Earth gravity 
model solutions from GRACE should especially improve the higher 
degree even zonal harmonics, so that it might happen that the differ- 
ence between the node-only LAGEOS-LAGEOS II and the node-only 
LAGEOS-LAGEOS II-Ajisai-Jason-1 combination will further enforce 
to the detriment of the former one, at least in regard to the gravita- 
tional error. The possibility of getting long time series of the Jason's 
node should seriously be investigated with real data tests. Moreover, 
Jason-1 is also affected by the orbital maneuvers but they are mainly 
in-plane. 
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